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ABSTRACT 

We  have  made  the  first  unambiguous  detection  of  vibrational  ground-state  maser  emission  from  ^*SiO  toward 
six  evolved  stars.  Using  the  Very  Large  Array,  we  simultaneously  observed  the  v  =  Q,J  =  1-0, 43.4-GHz,  ground- 
state  and  the  v  =  1,J=  1-0,  43. 1  -GHz,  first  excited-state  transitions  of  SiO  toward  the  oxygen-rich  evolved  stars 
IRC-i-1001 1,  o  Ceti,  W  Hya,  RX  Boo,  NML  Cyg,  and  R  Cas  and  the  S-type  star  x  Cyg.  We  detected  at  least  one 
V  =  0  SiO  maser  feature  from  six  of  the  seven  stars  observed,  with  peak  maser  brightness  temperatures  ranging 
from  10000  K  to  108800  K.  In  fact,  four  of  the  seven  v  =  0  spectra  show  multiple  maser  peaks,  a  phenomenon 
which  has  not  been  previously  observed.  Ground-state  thermal  emission  was  detected  for  one  of  the  stars,  RX 
Boo,  with  a  peak  brightness  temperature  of  200  K.  Comparing  the  v  =  0  and  the  v  =  1  transitions,  we  find  that  the 
ground-state  masers  are  much  weaker  with  spectral  characteristics  different  from  those  of  the  first  excited-state 
masers.  Eor  four  of  the  seven  stars  the  velocity  dispersion  is  smaller  for  the  v  =  0  emission  than  for  the  v  =  1 
emission,  for  one  star  the  dispersions  are  roughly  equivalent,  and  for  two  stars  (one  of  which  is  RX  Boo)  the 
velocity  spread  of  the  v  =  0  emission  is  larger.  In  most  cases,  the  peak  flux  density  in  the  v  =  0  emission  spectrum 
does  not  coincide  with  the  v  =  1  maser  peak.  Although  the  angular  resolution  of  these  VLA  observations  were 
insufficient  to  completely  resolve  the  spatial  structure  of  the  SiO  emission,  the  SiO  spot  maps  produced  from  the 
interferometric  image  cubes  suggest  that  the  v  =  0  masers  are  more  extended  than  their  v  =  1  counterparts. 

Subject  headings:  circumstellar  matter  —  masers  —  stars:  AGB  and  post-AGB 


1.  INTRODUCTION 

Although  the  study  of  SiO  maser  emission  toward  long- 
period  variable  stars  is  an  active  field  of  research,  very  little 
work  has  been  done  with  regards  to  the  vibrational  ground-state 
(v  =  0)  SiO  emission  toward  these  evolved  objects.  The  question 
of  the  existence  of  observable  ground-state  SiO  masers  was  first 
posed  by  Buhl  et  al.  (1975).  The  ground-state  lines  of  ^*SiO  and 
the  much  rarer  ^®SiO  and  ^°SiO  isotopic  species  have  been  ob¬ 
served  in  a  number  of  single-dish  surveys  conducted  primarily 
at  43  GHz  (v  =  0,7  =  1-0)  and  86  GHz  (v  =  0,7  =  2-  1).  In 
the  two  isotopic  species,  this  ground-state  emission  has  been 
shown  to  be  relatively  weak  unsaturated  maser  emission  (Al- 
colea  &  Bujarrabal  1992).  Eor  ^*SiO  the  nature  of  the  v  =  0 
emission  is  not  as  clear.  Although  Buhl  et  al.  (1975)  may  have 
detected  a  weak  v  =  0, 7  =  2-1  maser  line  toward  the  star  VY 
CMa,  the  general  consensus  on  the  86.8  GHz  line  is  that  it  is 
thermal  circumstellar  SiO  emission  with  broad  parabolic  or  flat- 
topped  line  profiles  (see  for  example  Bujarrabal,  et  al.  1986). 
The  v  =  0,7  =  1-0,  43.4-GHz  line  was  similarly  thought  to  be 
thermal  in  nature,  however  in  surveys  by  Jewell  et  al.  (1985) 
and  Jewell  et  al.  (1991)  several  of  the  v  =  0  line  profiles  showed 
narrow  emission  spikes  (<  1  km  s“')  superposed,  in  some  cases, 
on  broader  pedestals  of  emission.  Jewell  et  al.  suggested  that 
these  narrow  spikes  could  be  maser  emission  with  linewidths 
comparable  to  or  narrower  than  the  thermal  Doppler  width.  The 
data  of  Jewell  et  al.  (1991)  also  suggested  a  correlation  between 
the  intensities  of  the  v  =  0  and  v  =  1  lines  providing  additional 
evidence  for  ground-state  maser  emission  and  the  possibility 


that  the  gas  producing  the  two  types  of  emission  is  spatially 
coincident. 

To  date,  very  little  interferometric  imaging  of  the  v  =  0  emis¬ 
sion  toward  SiO  maser  emitting  regions  has  been  conducted. 
The  more  abundant  ^*SiO  species  was  observed  toward  the 
Orion  IRc2  star-forming  region  in  the  7=1-0  transition  (Chan¬ 
dler  &  DePree  1995)  and  the  7  =  2-1  transition  (Wright  et  al. 
1995).  In  each  case,  the  detected  emission  was  found  to  have 
both  maser  and  thermal  components.  Observations  of  the  rarer 
^®SiO  v  =  0, 7  =  2-1  transition  toward  Orion  (Baudry,  Herpin 
&  Lucas  1998)  also  showed  this  blend  of  thermal  and  maser 
emission.  A  comparison  of  the  ^*SiO  v  =  1, 7  =  2-1  and  ^®SiO 
v  =  0, 7  =  2-1  maser  spots  maps  lead  Baudry,  Herpin  &  Lucas 
(1998)  to  conclude  that  the  two  isotopic  species  are  not  excited 
in  the  same  gas  layers. 

Toward  late-type  stars,  interferometric  observations  of  the 
v  =  0, 7  =  2-1  ^*SiO  emission  have  been  conducted  by  Lucas 
et  al.  (1992)  and  by  Sahai  &  Bieging  (1993).  No  evidence  of 
ground-state  masers  was  found  in  either  observation  with  the 
possible  exception  being  an  intense  spike  of  emission  toward 
the  star  R  Cas  Lucas  et  al.  (1992).  Prior  to  the  work  presented 
here,  no  interferometric  studies  of  the  v  =  0,7  =  1  -  0  ^*SiO  tran¬ 
sition  had  been  conducted.  We  therefore  sought  to  detect  and 
image  the  ground-state  SiO  maser  emission  toward  7  stars  from 
the  list  of  possible  detections  in  Jewell  et  al.  (1991). 
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14.  ABSTRACT 

We  have  made  the  first  unambiguous  detection  of  vibrational  ground-state  maser  emission  from  28SiO 
toward  six  evolved  stars.  Using  the  Very  Large  Array,  we  simultaneously  observed  the  v  =  0,J  =  1-0, 
43.4-GHz,  groundstate  and  the  v  =  1,J  =  1-0,  43.1-GHz,  first  excited-state  transitions  of  28SiO  toward  the 
oxygen- rich  evolved  stars  IRC-f10011,  o  Ceti,  W  Hya,  RX  Boo,  NML  Cyg,  and  R  Cas  and  the  S-type  star  % 
Cyg.  We  detected  at  least  one  v  =  0  SiO  maser  feature  from  six  of  the  seven  stars  observed,  with  peak 
maser  brightness  temperatures  ranging  from  10000  K  to  108800  K.  In  fact,  four  of  the  seven  v  =  0  spectra 
show  multiple  maser  peaks,  a  phenomenon  which  has  not  been  previously  observed.  Ground-state  thermal 
emission  was  detected  for  one  of  the  stars,  RX  Boo,  with  a  peak  brightness  temperature  of  200  K. 
Comparing  the  v  =  0  and  the  v  =  1  transitions,  we  find  that  the  ground-state  masers  are  much  weaker  with 
spectral  characteristics  different  from  those  of  the  first  excited-state  masers.  Eor  four  of  the  seven  stars  the 
velocity  dispersion  is  smaller  for  the  v  =  0  emission  than  for  the  v  =  1  emission,  for  one  star  the  dispersions 
are  roughly  equivalent,  and  for  two  stars  (one  of  which  is  RX  Boo)  the  velocity  spread  of  the  v  =  0  emission 
is  larger.  In  most  cases,  the  peak  flux  density  in  the  v  =  0  emission  spectrum  does  not  coincide  with  the  v  = 

1  maser  peak.  Although  the  angular  resolution  of  these  VLA  observations  were  insufficient  to  completely 
resolve  the  spatial  structure  of  the  SiO  emission,  the  SiO  spot  maps  produced  from  the  interferometric 
image  cubes  suggest  that  the  v  =  0  masers  are  more  extended  than  their  v  =  1  counterparts. 
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Observations  of  the  v  =  0,J  =  1-0,  43.4-GHz  and  v  =  1,7  = 
1-0,  43.1-GHz  transitions  of  ^^SiO  were  performed  using  the 
Very  Large  Array  (VLA)  in  C-configuration.  The  VLA  is  main¬ 
tained  and  operated  by  the  National  Radio  Astronomy  Obser¬ 
vatory  (NRAO)*.  We  observed  the  7  stars  and  associated  extra- 
galactic  phase  calibrators  over  three  epochs.  The  stars  o  Ceti 
and  IRC H- 100 11  (WX  Psc)  were  observed  on  August  4,  2001; 
the  stars  NML  Cyg,  R  Cas,  and  x  Cyg  were  observed  on  Au¬ 
gust  8,  2001;  and  the  stars  RX  Boo  and  W  Hya  were  observed 
on  August  10,  2001.  Data  were  simultaneously  in  two  band- 
passes  of  bandwidth  6.25  MHz  centered  on  line  rest  frequencies 
of  43122.808  MHz  and  43423.858  MHz  for  the  v  =  1,7=  1-0 
and  v  =  0,7  =  1-0  transitions  respectively.  The  VLA  correlator 
in  normal  mode  produced  64  channels  per  IF  with  a  channel 
spacing  of  97.656  kHz  (^0.68  km  s“'). 

Data  were  reduced  using  the  standard  routines  within  the 
Astronomical  Image  Processing  System  (AIPS).  For  the  first 
epoch,  the  absolute  flux  density  scale  was  established  using 
3C84  assuming  a  flux  density  for  this  source  of  10.5  Jy.  For  the 
second  and  third  epochs,  both  3C84  and  3C286  were  used  for 
the  amplitude  calibration.  The  43-GHz  flux  density  of  3C286 
was  assumed  to  be  1 .49  Jy.  For  each  star  observed,  a  nearby 
extragalactic  calibrator  was  also  observed  in  order  to  estimate 
the  instrumental  and  atmospheric  phase  fluctuations.  Phase  cor¬ 
rections  estimated  from  the  calibration  sources  were  applied  to 
the  target  source  data.  To  eliminate  ringing  in  channels  adja¬ 
cent  to  strong  maser  features,  the  data  were  Hanning  smoothed 
prior  to  imaging.  This  resulted  in  a  degradation  of  the  spectral 
resolution  by  a  factor  of  2  to  (~  1.4  km  s“'). 

For  each  source,  an  iterative  self-calibration  and  imaging 
procedure  was  performed  to  map  a  single  strong  reference  spec¬ 
tral  feature  in  the  v  =  1  transition.  The  resulting  phase  solutions 
were  applied  to  all  channels  in  the  the  bands  containing  both  the 
V  =  1  and  V  =  0  transitions.  From  the  calibrated  data,  total  inten¬ 
sity  image  cubes,  consisting  of  256  x  256  pixel  images  of  the 
center  60  spectral  channels,  were  produced  for  both  the  v  =  1 
and  V  =  0  SiO  emission  lines  for  each  star.  Naturally- weighted 
individual  channel  images  were  made  with  an  approximately 
23"  X  23"  field  of  view.  The  synthesized  beam  sizes  ranged 
from  0.56"  x  0.54"  for  the  highest  declination  source,  R  Cas, 
to  1.09"  X  0.43"  for  the  lowest  declination  source,  W  Hya.  For 
the  star  RX  Boo,  the  v  =  0  emission  was  resolved  out  on  long 
baselines.  We  therefore  reduced  the  resolution  for  RX  Boo 
by  a  factor  of  two  by  limiting  the  uv  range  used  in  the  imag¬ 
ing.  This  resulted  in  an  image  cube  with  a  synthesized  beam  of 
1.12"  X  1.09" 

Positions  for  5  of  the  7  stars  as  determined  from  our  VLA 
observations  are  reported  in  Table  1 .  Each  position  was  deter¬ 
mined  from  a  two-dimensional  (2-D)  Gaussian  fit  to  the  peak  in 
the  image  of  the  reference  feature  for  the  v  =  1  transition  prior  to 
self-calibration.  The  errors  in  the  position  are  based  on  the  size 
of  the  synthesized  beam  divided  by  twice  the  signal-to-noise 
ratio  (SNR)  in  the  image.  These  errors  represent  a  small  frac¬ 
tion  of  the  beam  and  are  probably  not  representative  of  the  true 
error  in  the  derived  position  of  the  star.  There  may  be  signif¬ 
icant  blending  of  multiple  maser  features  in  the  peak  channel 
(discussed  below).  In  addition,  SiO  masers  typically  exhibit 
ring-like  structures  with  shell  sizes  are  on  the  order  of  tens  of 
milliarcseconds.  The  position  of  one  or  more  maser  features  at 
random  locations  on  the  ring  is  therefore  not  a  good  representa¬ 
tion  of  the  position  of  the  stellar  photocenter  and  may  be  off  by 


more  than  half  the  diameter  of  the  SiO  maser  shell.  For  the  stars 
RX  Boo  and  W  Hya,  the  phase  stability  prior  to  self-calibration 
was  insufficient  to  determine  a  reliable  position  from  the  image. 
Therefore,  no  position  is  reported  for  these  stars. 

Numerous  Very  Long  Baseline  Interferometry  (VLBI)  stud¬ 
ies  of  circumstellar  SiO  masers  have  shown  that  a  single  spec¬ 
tral  channel  often  contains  multiple  maser  features;  however, 
the  angular  resolution  of  our  VLA  data  was  not  sufficient  to  re¬ 
solve  multiple  features  separated  by  distances  less  than  typical 
maser  shell  diameters  of  30-50  mas.  In  order  to  parameterize 
individual  maser  components,  2-D  Gaussian  functions  were  fit 
to  the  peak  in  each  spectral  channel  of  the  image  cubes  using 
the  AIPS  task  SAD.  The  analysis  of  the  spectral  image  cubes, 
the  results  of  the  Gaussian  fits  to  the  v  =  0  and  v  =  1  SiO  emis¬ 
sion  and  the  astrophysical  interpretation  are  discussed  below. 


3.  RESULTS 


3.1.  Ground-State  Maser  Emission 


As  mentioned  in  the  previous  section,  we  fit  2-D  Gaussians 
to  the  emission  peak  in  each  spectral  channel  yielding  com¬ 
ponent  flux  density,  spatial  position,  and  velocity  information. 
Since  the  primary  goal  of  these  observations  was  to  deter¬ 
mine  whether  ground-state  SiO  maser  emission  is  exhibited  in 
the  atmospheres  of  late-type  stars,  we  concentrate  here  on  the 
ground-state,  v  =  0,  transition.  In  order  to  determine  whether 
the  observed  v  =  0  SiO  emission  has  a  maser  component  to  it, 
we  converted  the  fitted  component  flux  densities  to  equivalent 
brightness  temperatures  using: 


Tt  = 


SuM 

2kVt, 


:936- 


(1) 


where,  for  the  right-hand  estimate.  Si,  is  the  flux  density  in  Jy 
and  6s  is  the  geometric  mean  of  the  fitted  source  size  in  arcsec- 
onds.  Figures  1-7  plot  the  flux  densities  and  equivalent  bright¬ 
ness  temperatures  for  all  7  stars  observed.  In  these  figures,  the 
solid  line  represents  the  fitted  component  flux  densities  in  Jan- 
skys  (Jy),  and  the  points  represent  the  equivalent  component 
brightness  temperatures.  For  unresolved  maser  spots,  we  as¬ 
sumed  source  sizes  of  0.5  times  the  beam.  The  points  are  scaled 
by  the  logarithm  of  the  flux  density,  and  are  color  coded  accord¬ 
ing  to  velocity  bin.  The  flux  density,  deconvolved  source  size, 
and  brightness  temperature  for  the  strongest  spectral  feature  to¬ 
ward  each  source  are  tabulated  in  Table  2. 

The  kinetic  temperature  of  the  gas  in  the  circumstellar  enve¬ 
lope,  7):, can  be  estimated  from 

Tk  =  u^r°-^  (2) 

fi* 


where  T*  is  the  stellar  photospheric  temperature,  R*  is  the  stellar 
radius  and  R  is  the  radial  distance  in  the  circumstellar  envelope 
(CSE)  (Danger  &  Watson  1984;  Herpin  &  Baudry  2000).  As¬ 
suming  an  estimate  of  T*  w  2500  K,  the  kinetic  temperature  in 
the  region  of  the  v  =  1  SiO  masers  is  ~1500  K.  We  used  this  ki¬ 
netic  temperature  to  compare  to  the  brightness  temperatures  of 
the  SiO  ground-state  emission,  as  a  nominal  test  to  determine 
whether  the  emission  observed  in  the  v  =  0,7  =  1-0  transition 
is  indeed  masing.  This  temperature  is  rather  conservative  given 
the  expectation  that  the  v  =  0  masers  should  occur  at  a  larger  ra¬ 
dial  distance  than  the  v  =  1  SiO  masers  (discussed  later).  Erom 
Table  2  we  see  that  for  6  of  the  7  stars  observed,  the  peak  bright¬ 
ness  temperature  is  much  higher  than  the  nominal  1500  K  ki¬ 
netic  temperature.  We  take  this  comparison  to  clearly  indicate 
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that  some  of  the  v  =  0,7  =  1-0  SiO  emission  is  indeed  maser 
emission.  The  only  star  that  was  observed  that  does  not  show 
maser  emission  from  ground-state  SiO  in  is  RX  Boo.  In  ad¬ 
dition,  the  deconvolved  angular  sizes  for  the  peaks  in  the  six 
stars  with  suspected  maser  emission  (see  Table  2)  were  well 
below  the  corresponding  beam  size.  This  is  also  suggestive 
of  compact  maser  emission.  The  only  source  with  a  decon¬ 
volved  source  size  greater  than  the  synthesized  beam  was  again 
RX  Boo. 

Examination  of  Eigures  1-7  shows  that  4  of  the  6  ground- 
state  maser  sources  show  multiple  brightness  temperature  peaks 
well  above  the  nominal  kinetic  temperature.  The  two  excep¬ 
tions  are  o  Ceti  and  x  Cyg.  It  is  highly  probably  that  these 
additional  features  are  also  ground-state  masers.  Such  multi¬ 
ple  masing  features  contradict  previously  reported  detections 
which  showed  single  narrow  emission  peaks  in  the  spectra.  The 
six  V  =  0  maser  stars  appear  to  show  a  range  in  the  degree  of 
masing  from  x  Cyg  and  o  Ceti  which  exhibit  single  maser  peaks 
to  sources  like  NML  Cyg  and  IRCh-1001  1  in  which  nearly  the 
entire  spectrum  is  above  the  brightness  temperature  threshold 
for  masing.  The  spectra  for  the  latter  two  stars  are  similar  to  the 
v  =  Q,J  =  1-0  spectrum  for  Orion  IRc2,  which  Chandler  &  De- 
Pree  (1995)  found  to  be  dominated  by  SiO  maser  emission.  Eor 
X  Cyg,  it  appears  that  the  single  maser  feature  sits  on  a  broader 
pedestal  of  thermal  SiO  emission.  This  is  reminiscent  of  some 
of  the  V  =  Q,J  =1-0  spectra  observed  by  Jewell  et  al.  (1991) 
and  the  v  =  0, 7  =  2-1  spectrum  for  R  Cas  observed  by  Lucas 
etal.  (1992). 

3.2.  Comparison  of  the  v  =  0  and  v  =  1  SiO  Emission 

Having  established  that  there  is  at  least  some  degree  of  mas¬ 
ing  in  the  ground-state  SiO  emission  toward  6  of  the  7  stars 
observed,  we  sought  to  compare  the  v  =  0  emission  with  that 
observed  for  the  more  common  v=l, 7=1-0,  first  excited-state 
maser  transition.  To  facilitate  this  comparison,  we  plotted  the 
results  of  the  2-D  Gaussian  fits  in  Eigures  8-14.  In  each  figure, 
the  two  left  panels  show  the  spectrum  (upper  panel)  and  spatial 
distribution  (lower  panel)  of  the  v  =  0  SiO  emission.  Similarly, 
the  two  right  panels  show  the  spectrum  and  spatial  distribution 
of  the  V  =  1  SiO  emission.  Components  derived  from  the  2-D 
fits  are  represented  by  color-coded  circles  with  sizes  scaled  by 
the  logarithm  of  the  flux  density  and  colors  corresponding  to 
particular  velocity  bins.  The  total  velocity  range  and  velocity 
binning  is  the  same  for  each  transition  for  each  star.  The  flux 
density  scale  in  the  upper  panels  and  the  right  ascension  and 
declination  offsets  in  the  lower  panels  are  adjusted  to  provide 
the  best  plot  filling  factor  in  each  transition.  Also  plotted  in  the 
upper  panels  of  each  figure,  are  solid  lines  which  represent  the 
vector-averaged,  cross-power,  spectrum  averaged  over  all  VLA 
antennas  for  comparison. 

Comparing  the  upper  panels  in  Eigures  8-14  we  notice  that 
the  spectra  for  the  two  SiO  transitions  toward  each  star  are  quite 
different.  Aside  from  the  obvious  flux  density  differences,  the 
spectra  vary  in  number  of  peaks,  center  velocities  of  the  peaks, 
and  total  velocity  dispersion  in  the  transition.  In  Table  3  we  list 
some  of  the  relevant  spectral  parameters  derived  from  the  fits 
to  the  image  data.  Both  the  figures  and  the  table  show  that  for 
majority  of  the  stars,  x  Cyg  being  the  one  exception,  the  peak 
flux  densities  in  the  two  transitions  do  not  coincide  in  velocity 
to  within  1  km  s“' .  Additionally,  the  velocity  dispersions  (AV 
in  Table  3)  between  the  two  transitions  vary  for  each  star.  In 
4  of  the  7  stars  AV  is  larger  for  the  v  =  1  transition  than  for 


the  V  =  0  line.  Two  of  the  stars,  IRC-tlOOl  1  and  RX  Boo,  have 
larger  dispersions  in  the  v  =  0  transition,  although  RX  Boo  is  not 
masing  in  this  transition,  and  one  star,  NML  Cyg,  has  velocity 
dispersions  which  are  roughly  equivalent. 

Plotted  in  the  lower  panels  of  Eigures  8-14  are  the  spatial  dis¬ 
tributions  of  the  SiO  emission  in  each  transition  for  each  star. 
As  mentioned  in  Section  ??,  the  VLA  in  C-configuration  is  in¬ 
sufficient  to  resolve  multiple  maser  components  within  a  spec¬ 
tral  channel,  thus  we  fit  a  single  2-D  Gaussian  to  the  emission 
peak  within  each  channel.  Eor  the  v  =  1  maser  emission,  the  po¬ 
sition  errors  reported  by  the  AIPS  task  SAD  for  these  fits  range 
from  ~15  mas  to  less  than  1  mas.  Eor  the  weaker  v  =  0  maser 
emission,  position  errors  ranged  from  ~70  mas  to  1  mas.  Eor 
the  thermal  v  =  0  emission  toward  RX  Boo,  errors  ranged  from 
100  mas  to  30  mas.  These  reported  errors  are  based  on  the  size 
of  the  beam  divided  by  twice  the  signal-to-noise  ratio  (SNR) 
in  the  channel  image  and  are  considered  more  realistic  than  the 
formal  errors  of  the  Gaussian  fit.  Nevertheless,  the  reported 
errors  are  still  a  very  small  fraction  (0.1-0.002)  of  the  synthe¬ 
sized  beam  and  some  may  question  their  reliability.  Because  of 
this  and  the  fact  that  there  is  probably  some  component  blend¬ 
ing  occurring  in  the  images,  the  spatial  distributions  mapped  in 
Eigures  8-14  are  discussed  only  in  a  qualitative  sense.  A  true 
quantitative  comparison  of  the  spatial  distributions  of  the  two 
masing  transitions  will  have  to  wait  for  higher-resolution  inter¬ 
ferometric  observations. 

With  this  in  mind,  we  note  that  for  all  of  the  stars  observed, 
the  spatial  extent  of  the  v  =  0  SiO  emission  is  several  times 
larger  than  that  of  the  corresponding  v  =  1  maser  emission.  In 
addition,  we  find  that  the  v  =  1  masers  form  compact  distri¬ 
butions  ranging  from  ~30-90  mas  consistent  with  typical  SiO 
maser  shell  sizes  as  measured  by  VLBI.  Especially  enticing  are 
the  maps  for  NML  Cyg  (see  Eigure  1 3)  in  which  we  find  that  the 
V  =  1  maser  distribution  agrees  with  previous  VLBI  results  from 
Boboltz  &  Marvel  (2000).  Additionally,  the  v  =  0  ground-state 
masers  appear  to  be  resolved  into  a  definite  ring-like  structure 
with  a  radius  of  about  50  mas.  Using  the  Berkeley  Infrared  Spa¬ 
tial  Interferometer  (ISI),  Danchi  et  al.  (2001)  determined  the 
radius  of  the  inner  dust  shell  of  NML  Cyg  to  be  80-100  mas. 
This  would  place  the  v  =  0  masers  in  a  region  roughly  half  the 
radial  distance  between  the  v  =  1,7  =  1-0  masers  and  the  inner 
circumstellar  dust  shell  for  NML  Cyg.  Assuming  a  stellar  size 
of  approximately  9  mas  (Boboltz  &  Marvel  2000),  the  kinetic 
temperature  of  the  SiO  in  the  v  =  0  region  is  ~750  K  based  on 
Equation  2. 

4.  DISCUSSION 

Ground-State,  ^^SiO  maser  emission  is  rarely  addressed  in 
the  literature  in  which  numerical  models  for  the  generation  of 
SiO  masers  are  presented.  This  may  be  due  to  the  fact  that 
to  date  there  has  been  very  little  evidence  for  the  existence  of 
masers  in  this  transition.  In  an  early  paper,  Kwan  &  Scov- 
ille  (1974)  commented  that  population  inversions  in  the  v  =  0 
states  are  possible  under  certain  conditions.  They  stated  that 
the  percentage  of  the  inversion  is  small  due  to  the  fact  that  the 
collision  rate  among  the  v  =  0  levels  is  comparable  to  the  rate 
of  radiative  depopulation.  In  one  model  consisting  of  strictly 
radiative  excitation,Kwan  &  Scoville  (1974)  found  an  inver¬ 
sion  of  the  V  =  Q,J  =  1-0  line  which  occurs  at  a  distance  of 
5  X  10'‘*-4  X  10^^  cm  (~9-69  Rf)  with  an  emission  rate  down 
by  a  factor  of  100  from  the  v  =  1,7  =  1  -0  line. 

More  recent  investigations  addressing  the  excitation  of  the 
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SiO  masers  through  overlaps  with  optical  (Rausch  et  al.  1996) 
and  infrared  (Herpin  &  Baudry  2000)  lines,  demonstrated  pos¬ 
sible  mechanisms  for  the  generation  of  v  =  0,  ^*SiO  masers. 
In  the  case  of  optical  pumping,  Rausch  et  al.  (1996)  found  the 
ground-state  emission  to  be  thermal  for  the  v  =  0, 7  =  2-1  and 
v-0,7  =  3-2  transitions  of  ^*SiO,  but  that  maser  emission 
was  generated  for  the  v  =  0,7  =  1-0  state  in  their  simulations. 
The  infrared  pumping  mechanism  of  Herpin  &  Baudry  (2000) 
was  used  to  resolve  some  of  the  previously  unexplained  detec¬ 
tions  of  certain  SiO  maser  transitions  and  to  predict  some  pos¬ 
sible  new  SiO  masers.  In  addition,  their  numerical  model  also 
demonstrated  that  v  =  Q,J  =  1-0  ^^SiO  could  be  generated.  In 
both  the  optical  and  infrared  pumping  cases,  the  v  =  0  masers 
are  generated  in  regions  of  the  CSE  where  the  SiO  abundances 
are  lower  than  for  the  typical  excited-state  masers.  Since  the 
SiO  abundance  decreases  with  distance  from  the  star  as  SiO 
condenses  onto  dust  grains,  then  the  v  =  0  masers  are  natu¬ 
rally  formed  at  distances  beyond  that  of  the  vibrationally  ex¬ 
cited  masers  in  these  models. 

Our  qualitative  comparison  of  the  spatial  distributions  sup¬ 
ports  the  idea  that  the  v  =  0  SiO  masers  occur  at  distances 
greater  than  v  =  1  maser  shell.  In  this  case  the  v  =  0  maser 
distribution  is  several  times  the  diameter  of  the  corresponding 
V  =  1  masers.  This  contradicts  the  speculation  by  Jewell  et  al. 
(1991)  that  the  correlation  between  the  intensities  of  the  v  =  0 
and  V  =  1  lines  suggests  that  emission  from  the  two  transitions 
is  spatially  coincident.  Baudry,  Herpin  &  Lucas  (1998),  in  their 
interferometric  observations  of  the  v  =  0,7  =  2- 1  ^®SiO  and  the 
v=l, 7  =  2-1  ^^SiO  masers  toward  Orion  IRc2,  determined  that 
the  ground-state  and  excited-state  masers  were  not  co-spatial. 
They  state  that  this  is  a  consequence  of  all  collisional/radiative 
pumping  models  which  predict  that  various  vibrational  state 
masers  should  peak  in  different  spatial  regions  of  the  CSE. 


5.  SUMMARY 

Interferometric  observations  using  the  VLA  have  allowed  us 
to  definitively  detect  ground-state  v  =  0,7  =  1-0  ^*SiO  masers 
toward  six  long-period  variable  stars,  IRSh-10011  (WX  Psc), 
o  Ceti,  W  Hya,  NML  Cyg,  R  Cas,  and  x  Cyg.  Thermal  v  =  0 
emission  was  detected  for  the  only  other  the  star  in  our  sam¬ 
ple,  RX  Boo.  Peak  brightness  temperatures  for  the  ground-state 
masers  ranged  from  10000  K  to  108800  K  at  the  time  of  our  ob¬ 
servations.  Of  the  six  ground-state  maser  stars,  four  exhibited 
multiple  maser  peaks,  a  phenomenon  not  previously  observed. 

Comparing  the  v  =  0  transition  with  the  simultaneously  ob¬ 
served  v  =  1,7  =  1-0  maser  emission,  we  find  that  the  ground- 
state  masers  are  weaker  with  spectra  different  than  the  corre¬ 
sponding  first  excited-state  masers.  In  general  the  peak  in  the 
V  =  0  spectrum  for  a  given  star  occurs  at  a  different  velocity 
from  that  of  the  v  =  1  maser  peak.  Additionally,  the  dispersions 
in  velocity  are  not  the  same. 

We  were  able  to  make  low-resolution  image  cubes  of  both 
transitions  for  each  star  from  our  VLA  C-configuration  data. 
Erom  these  images  we  created  spot  maps  of  the  emission.  Al¬ 
though  the  resolution  of  our  observations  was  insufficient  to 
provide  a  quantitative  comparison  of  the  spatial  morphologies 
of  the  two  transitions,  we  can  say  that  the  ground-state  SiO 
masers  appear  to  be  distributed  over  a  larger  region  of  the  CSE 
than  their  v  =  1  counterparts.  Euture  higher-resolution  obser¬ 
vations  should  enable  us  to  clearly  resolve  the  spatial  distribu¬ 
tion  of  ground-state  masers  toward  these  long-period  variables. 
The  precise  placement  of  these  masers  in  the  circumstellar  en¬ 
velope  relative  to  the  v  =  1  masers  and  to  the  circumstellar  dust 
shell  has  important  implications  for  circumstellar  chemistry  and 
maser  pumping  mechanisms. 
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Table  1 

Positions  of  v  =  1,7  =  1  -0  SiO  maser  reference  eeatures. 


Source 

^peak. 

(km  ) 

Q  (J2000) 

(h  m  s) 

&  (J2000) 

(0  , 

IRC+IOOII 

10.7 

01  06  25.9886  ±0.00004  (±0.0006") 

12  35  53.0373  ±0.0006 

o  Ceti . 

44.3 

02  19  20.7909  ±0.00003  (±0.0004") 

-02  58  39.8577  ±0.0005 

W  Hya  . . . . 

40.7 

RX  Boo  . . . 

0.3 

xCyg  .... 

9.7 

19  50  33.9177  ±0.00022  (±0.0028") 

32  54  50.5987  ±0.0021 

NMLCyg  . 

-15.2 

20  46  25.5430  ±0.00007  (±0.0008") 

40  06  59.4177  ±0.0010 

R  Cas . 

25.4 

23  58  24.8805  ±0.00015  (±0.0014") 

51  23  19.7645  ±0.0014 

Table  2 

Peak  elux  densities  and  brightness  temperatures  oe  v  =  0,7  =  1  -  0  SiO  emission. 


Source 

Peak  Flux 
Density 

(Jy) 

Velocity 
(km  s“‘) 

Angular 

Size 

(") 

Brightness 

Temp. 

(K) 

IRC+lOOll 

2.90±  0.02 

8.0 

0.16 

108800 

0  Ceti . 

1.68±  0.02 

46.4 

0.26 

23400 

W  Hya  ... . 

3.74±  0.05 

36.0 

0.42 

19900 

RX  Boo  . . . 

1.81±  0.05 

3.7 

2.34 

200 

xCyg  .... 

1.24±  0.02 

9.0 

0.34 

10000 

NML  Cyg  . 

4.18±0.03 

-1.0 

0.19 

105900 

R  Cas . 

4.54±  0.04 

30.8 

0.20 

107600 

Table  3 

Observed  spectral  properties  of  the  v  =  0  and  v  =  1  SiO  emission. 


Source 

AV 

(km  ) 

v  =  0,7  =  1 

kpeak 
(km  S“^) 

-0 

Peak  Flux 
Density 

(Jy) 

AV 

(km  ) 

v=  1,7=  1 

kpeak 
(km  S“^) 

-0 

Peak  Flux 
Density 

(Jy) 

IRC+lOOll 

15.5 

8.0 

2.90±  0.02 

12.2 

10.7 

298.0  ±  0.6 

o  Ceti . 

7.4 

46.4 

1.68±  0.02 

16.3 

44.3 

617.6  ±0.9 

W  Hya  ... . 

11.5 

36.0 

3.74±  0.05 

17.7 

40.7 

551.1  ±4.4 

RX  Boo  . . . 

16.2 

3.7 

1.1 8±  0.05 

11.5 

0.3 

20.1  ±0.2 

xCyg  .... 

12.8 

9.0 

1.24±  0.02 

15.6 

9.7 

19.25±  0.05 

NML  Cyg  . 

35.1 

-1.0 

4.18±  0.03 

36.0 

-15.2 

12.25±  0.04 

R  Cas . 

12.8 

30.8 

4.54±  0.04 

17.0 

25.4 

475.9  ±  1.9 
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Fig.  8. —  Spatial/velocity  structure  of  the  v  =  0,7  =  1-0  (left  panels)  and  v  =  1,7  =  1  -0  (right  panels)  SiO  emission  toward  IRC+1001 1  (WX  Psc)  as  measured 
by  the  2-D  Gaussian  fits  to  the  peaks  in  the  VLA  channel  images.  The  two  upper  panels  show  the  spectra  formed  by  plotting  flux  density  versus  LSR  velocity,  color 
coded  in  velocity  increments  from  blue-shifted  (left)  to  red-shifted  (right).  The  two  lower  panels  are  spatial  spot  maps  of  the  of  the  SiO  emission.  The  color  of  each 
point  is  representative  of  the  con'esponding  velocity  bin  in  the  spectrum  and  the  size  of  each  point  is  proportional  to  the  logarithm  of  the  flux  density.  The  solid  lines 
in  the  upper  panels  represent  the  vector-averaged,  cross-con'elated  flux  density  in  each  channel  averaged  over  all  antennas. 
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Fig.  9. —  Same  as  Figure  8  for  star  o  Ceti. 
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Fig.  10. —  Same  as  Figure  8  for  star  W  Hya. 
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Fig.  1 1 . —  Same  as  Figure  8  for  star  RX  Boo. 
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Fig.  12. —  Same  as  Figure  8  for  star  x  Cyg. 
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Fig.  13. —  Same  as  Figure  8  for  star  NML  Cyg. 
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Fig.  14. —  Same  as  Figure  8  for  star  R  Cas. 


